Dynamics of DNA-breathing: Weak noise analysis, finite time singularity, and 
mapping onto the quantum Coulomb problem 
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We study the dynamics of denaturation bubbles in double-stranded DNA on the basis of the 
Poland- Scheraga model. We show that long time distributions for the survival of DNA bubbles and 
the size autocorrelation function can be derived from an asymptotic weak noise approach. In partic- 
ular, below the melting temperature the bubble closure corresponds to a noisy finite time singularity. 
We demonstrate that the associated Fokker-Planck equation is equivalent to a quantum Coulomb 
problem. Below the melting temperature the bubble lifetime is associated with the continuum of 
scattering states of the repulsive Coulomb potential; at the melting temperature the Coulomb po- 
tential vanishes and the underlying first exit dynamics exhibits a long time power law tail; above 
the melting temperature, corresponding to an attractive Coulomb potential, the long time dynamics 
is controlled by the lowest bound state. Correlations and finite size effects are discussed. 

PACS numbers: 05.40.-a,02.50.-r,87.15.-v,87.10.-|-e 



I. INTRODUCTION 

Under physiological conditions the Watson-Crick 
double-helix of DNA constitutes the equilibrium struc- 
ture, its stability ensured by hydrogen-bonding of paired 
bases and base stacking between nearest neighbor pairs of 
base pairs [H, |1] ■ By variation of temperature or pH- value 
double-stranded DNA progressively denatures, yielding 
regions of single-stranded DNA, until the double-strand 
is fully molten. This is the helix-coil transition taking 
place at a melting temperature Tm defined as the temper- 
ature at which half of the DNA molecule has undergone 
denaturation 

However, already at room temperature thermal fluc- 
tuations cause rare opening events of small denaturation 
zones in the double-helix [J • These DNA bubbles consist 
of flexible single-stranded DNA, and their size fluctuates 
in size by step-wise zipping and unzipping of the base 
pairs at the two zipper forks where the bubble connects 
to the intact double-strand. Below the melting tempera- 
ture Tm, once formed, a bubble is an intermittent feature 
and will eventually zip close again. The multistate DNA 
breathing can be monitored in real time on the single 
DNA level 0]. Biologically, the existence of intermittent 
(though infrequent) bubble domains is important, as the 
opening of the Watson-Crick base pairs by breaking of 
the hydrogen bonds between complementary bases dis- 
rupts the helical stack. The flipping out of the ordered 
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stack of the unpaired bases allows the binding of specific 
chemicals or proteins, that otherwise would not be able 
to access the reactive sites of the bases H, S B Q • 

The size of the bubble domains varies from a few bro- 
ken base pairs well below Tm, up to some two hundred 
closer to Tm- Above T,„, individual bubbles continuously 
increase in size, and merge with vicinal bubbles, until 
complete denaturation Assuming that the bubble 
breathing dynamics takes place on a slower time scale 
than the equilibration of the DNA single-strand consti- 
tuting the bubbles, DNA-breathing can be interpreted 
as a random walk in the ID coordinate x, the number of 
denatured base pairs. 

DNA breathing has been investigated in the Dauxois- 
Peyrard-Bishop model [1, Q , that describes the motion 
of coupled oscillators representing the base pairs. On 
the basis of the Poland-Scheraga model, DNA breathing 
has been studied in terms of continuous Fokker-Planck 
approaches [13, [Hi , and in terms of the discrete master 
equation and the stochastic Gillespie scheme [H, [H, 0, 
[isl . [igI . [T3| ■ The coalescence of two bubble domains was 
analyzed in Ref. p^ . 

In what follows we study the Langcvin and Fokker- 
Planck non-equilibrium extension of the Poland-Schcraga 
model in terms of both a general weak noise approach ac- 
cessing the long time behavior, see e.g. , Refs. Il9l . [20i , and 
a mapping to a quantum Coulomb problem j2l|. This al- 
lows us to investigate in more detail the finite time singu- 
larity underlying the breathing dynamics, as well as the 
survival of individual bubbles. The paper is organized in 
the following manner. In Sec.|TTl we introduce and discuss 
the model, in Sec. IIIII we apply the weak noise approach 
and extract long time results and study the stability of 
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the solutions. In Sec. lIVI wc map the problem to a quan- 
tum Coulomb problem and derive the long-time scaling 
of the bubble survival. Finally, in Sec. |V] we discuss the 
results and draw our conclusions in Sec. I VII 

II. DYNAMIC MODEL FOR DNA BREATHING 

In the Poland-Scheraga free energy approach, bubbles 
are introduced as free energy changes to the double- 
helical ground state, such that the disruption of each 
additional base pair of a bubble requires to cross an ener- 
getic barrier that is rewarded by an entropy gain. While 
the persistence length of double-stranded DNA is rather 
large (of the order of 50nm) and it is assumed to have no 
configurational entropy, the single-stranded bubbles are 
flexible, and therefore behave like a polymer ring. The 
Poland-Scheraga partition factor for a single bubble in a 
homopolymer is of the form 

Z(m) =ao""(l + m)-", (2.1) 

where m counts the (discrete) number of broken base 
pairs, and u = exp (— /37), with (3 = l/[fcT], is the Boltz- 
mann factor for breaking the stacking interactions when 
disrupting an additional base pair. The cooperativity 
factor (To = exp (— /37o) quantifies the so-called boundary 
energy 70 for initiating a bubble. 70 is of the order of 
8000 cal/mol, corresponding to approximately 13 kT at 
37° C. Occasionally, somewhat smaller values for are 
assumed, down to approximately 8 kT . Bubbles below 
the melting point of DNA are therefore rare events. Typ- 
ical equilibrium melting temperatures of DNA for stan- 
dard salt conditions are in the range ^ 70 — 100°C, 
depending on the relative content of weaker AT and 
stronger GC Watson-Crick base pairs. Thus, double- 
stranded DNA denatures at much higher temperatures 
as many proteins. Note that the melting temperature 
of DNA can also be increased by change of the natural 
winding, as opening of the double-strand in ring DNA is 
coupled with the creation of superstructure; this is the 
case, for instance, in underwater bacteria living in hot 
vents, compare Ref. [2^ , and references therein. 

Due to the large value of (Tq, below the melting tem- 
perature to good approximation individual bubbles are 
statistically independent, and therefore a one-bubble pic- 
ture appropriate. Having experimental setups in mind 
as realized in Ref. where special DNA constructs are 
designed such that they have only one potential bubble 
domain, we also consider a one-bubble picture at and 
above T^. Our results are meant to apply to such typ- 
ical single molecule setups. In comparison to the rather 
high energy barrier 70, according to which the opening 
of a bubble corresponds to a nucleation process, to break 
the stacking of a single pair of base pairs requires much 
less thermal activation, ranging from 7 = —0.1 to -1-3.9 
kT for TA/AT and GC/CG pairs of base pairs at 37°C, 
respectively; here, the positive sign refers to a thermody- 
namically stable state. These comparatively low values 



for the stacking free energy of base pairs stems from the 
fact that stacking enthalpy cost and entropy release on 
base pair disruption almost cancel. Finally, the term 
(1 -|- m)~'^ measures the entropy loss on formation of a 
closed polymer ring, with respect to a linear chain of 
equal length. The offset by 1 is often taken into ac- 
count to represent the short persistence length of sin- 
gle stranded DNA. For the critical exponent c, one typ- 
ically uses the value 1.76 of a Flory chain in three di- 
mensions 0, [H, m, [13, m, [2^, while a slightly larger 
value (c = 2.12) was suggested based on different poly- 
mer models [13, [H, ^^M, IM [H]. Here, we disregard 
the offset, and consider the pure power-law form m"*^. 

In the following, we consider the continuum limit of 
the above picture, measuring the "number" of broken 
base pairs with the continuous variable x. The Poland- 
Scheraga free energy for a single bubble then has the form 

= -fn+-fx + ckT\nx. (2.2) 

where cc > is the bubble size as measured in units 
of base pairs. Treating the bubble size x as a contin- 
uum variable, we impose an absorbing wall at a; = 0, the 
zero-size bubble. The completely closed bubble state is 
stabilized by the size of the cooperativity factor ctq > £^nd 
bubbles therefore become rare events. Expression (|2.2p 
corresponds to a logarithmic sink in ^ at x = 0. The 
free energy density 7(T) has a temperature dependence, 
which we write as 

7(r) = 7i(T„,-r)/r™, (2.3) 

where T,„ is the melting temperature. 

From Eq. (|2.2p it follows that a characteristic bubble 
size is set by xi = ckT/\'y\. For large bubble size x > xi 
the linear term dominates and the free energy grows like 
=^ ~ 7o + 72^- For small bubbles x < xi [or close to T™, 
where 7(r) ~ 0] the free energy is characterized by the 
logarithmic sink but has strictly speaking a minimum 
at ^ = 7o for zero bubble size. We distinguish two 
temperature ranges: 

(i) For 7 < 0, i.e., T > Tm, the free energy has a maxi- 
mum ,^max = 7o + cA:T(loga;i — 1) aX x = xi. The free en- 
ergy profile thus defines a Kramers escape problem in the 
sense that an initial bubble can grow in size correspond- 
ing to the complete denaturation of the double stranded 
DNA. The escape probability Pose oc exp(-A^/fcr), 
where the free energy barrier is = ckT{\ogxi — 1), 
i.e., 

Pe.c«(^)"^ (2.4) 

(ii) For 7 > 0, i.e., T < T,n, the free energy increases 
monotonically from = 70 at x = and the finite size 
bubbles are stable. The change of sign of 7 at T = T^ 
thus defines the bubble melting. 

For 7 < 0, i.e., T > Tm, the free energy has a maximum 
and decreases for large bubble size, as a result the bubbles 
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energy (|2.2p in Eq. (|2.5p we have in particular 
at X 



(2.7) 



where we have found it convenient to introduce the in- 
verse time scales J^i and Vi2, 

VLi = Dc/sT, (2.8a) 
n2 = -D^ = D-ii{T-T„,)/T„,. (2.8b) 

Note that the characteristic bubble size xi = ckT/\^\ is 
given by 



Xi 



ckT 

1 



(2.9) 



and thus emerges from the time scale competition be- 
tween the VLi , from a dynamic point of view. 

In the limits of large and small bubble sizes, the 
Langevin equation (|2.5[) allows exact solutions: 

(i) For large bubble size x 3> xi we can ignore the loop 
closure or entropic contribution ckT/x and wo obtain the 
Langevin equation 



dt 



(2.10) 



describing a ID random walk with an overall drift ve- 
locity Q2- For large x we thus obtain the distribution 



Figure 1: We depict the free energy profile =^ — 70 below and 
above the melting temperature Tm as a function of bubble 
size. In a) we show =^ — 70 for 7 > 0, i.e., T < Tm; in b) 
we show — 70 for 7 < 0, i.e., T > Tm- For large bubble 
sizes, X ^ xi the free energy behaves approximately linearly 
as function of bubble size. For small bubble sizes the free 
energy has a logarithmic sink corresponding to the absorbing 
state at X = (arbitrary units). Above melting, there exists a 
nucleation barrier that needs to be crossed before the bubble is 
allowed to grow towards full denaturation. In both cases, the 
comparatively high initiation barrier 70 has to be overcome 
to seed the bubble. 



expand and the double stranded DNA denatures, that is, 
melts. In Fig. [T] we have depicted the free energy profile 
as a function of bubble size for 7 > 0, T < Tm, and for 
7 < 0, T > Tm. 

The stochastic bubble dynamics in the free energy 
landscape ^ is described by the Langevin equation 



dx _ ^d^ 

dt dx ' 



driven by thermal noise ^, that is characterized by the 
correlation function 



m^{t'))^'iDkT5{t~t'). 



(2.6) 



The kinetic coefficient D of dimension [kT) ^ sets 
the inverse time scale of the dynamics. Inserting the free 



P{x,t) 



1 



ViirDkTt 



exp 



(x - xo - ii.2ty 



ADkTt 



(2.11) 



where xq is the initial (large) bubble size. It follows that 
the mean bubble size scales linearly with time, (x) ~ xo + 
il2t. Below T,n {^2 < 0) the bubble size shrinks towards 
bubble closure; above T,n {^2 > 0) the bubble size grows, 
leading to denaturation. The mean square bubble size 
fiuctuations ((Ax)^) = 2DkTt, increase linearly in time, 
a typical characteristic of a random walk. 

Taking into account the absorbing state condition 
P(x = 0,t) ~ for zero bubble size by forming the lin- 
ear combination (method of images), we obtain for the 
distribution [3J| 



P. 



1 



abs 



exp 



exp 



^/AnDkTt 
X0O2 1 



(x - xq - n2tY 

WkTt 

(x -I- Xo — ^2t)^ 



ADkTt 



(2.5) and infer, using the definition [si 



W{t)^ 

the first passage time density 



, 9Pabs 

dt 



W{t) 



Xo 



y/AirDkTt^ 



exp 



(xq + Vl2tf 

ADkTt 



(2.12) 



(2.13) 



(2.14) 
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Figure 2: In a) we show the time evolution of a small bubble of 
size X in the absence of thermal noise. For a:: = correspond- 
ing to bubble closure we encounter a finite-time-singularity 
at to = xo/2Q,\. In b) we depict the noisy case. Here the 
first passage time is a statistical event characterized by W{t) 
(arbitrary units). 



with the typical Sparre Andersen asymptotics 



W{t) 



zt 



y/4:TTDkT 



-3/2 



(2.15) 



(ii) For small bubble size a; <C xi the nonlinear entropic 
term dominates and the bubble dynamics is governed by 
the nonlinear Langevin equation 



P(0,t) = 0, has the form 

^ni/2DkT+l/2 ^-(x^+xl)/4DkTt 



P{x,t) 



X, 



ni/DkT-l/2 







/ xxo 

X/l/2+^2,/2DfcT^2^^y 



2DkTt 

V 



(2.17) 



Here Ii, is the Bessel function of imaginary argument, 
Iv{z) = {—lY Jv{iz) d^. Correspondingly, we find the 
first passage time distribution 



W{i) = 



■ exp 



r(l/2-Oi/2L>A:T) 
x(4i?m)-3/2-f^i/2i^^T 

with the long time tail 



ADkTt J 



(2.18) 



Wit) 



l+Oi/DfeT 3/2-C/2 



r(l/2 - 17i/2L>fcr)(4L>fcr)i/2+ni/2DfeT ' 

(2.19) 

where we substituted back for ili: For small bubble sizes, 
the exponent c due to the polymeric interactions changes 
the first passage statistics. As already noted in Ref. [2l| . 
this modified exponent for c > 1 gives rise to a finite 
mean first passage time tW{t)dt, in contrast to the 
first passage time distribution (|2.14|) 

In the general case for bubbles of all sizes the fluctua- 
tions of double-stranded DNA is described by Eq. (|2.7p . 
The associated Fokkcr-Planck equation for the distribu- 
tion P{x, t) has the form (compare also Refs. [O, HT], [2^ ) 



dP 
'dt 



d_ 

dx 



X 



P + DkT 



d^P 



dx"^ 



(2.20) 



and provides the complete description of the single bub- 
ble dynamics in double-stranded homopolymer DNA in 
the continuum limit of the Poland-Scheraga model. For 
large bubble sizes where the entropic term f^i/a; can be 
neglected the solution of Eq. (|2.20p is given by Eqs. (|2.1ip 
and (|2.12p . Conversely, for small bubble sizes, where the 
entropic term Q,i/x dominates, or for all bubble sizes pre- 
cisely at the transition temperature f22 = (T = Tm), 
the solution of Eq. (|2.20p is given by the noisy flnite-time- 
singularity solution in Eqs. (|2.17[) and (|2.18p . 



dx 



(2.16) 



For vanishing noise Eq. (|2.16p has the solution x = 
(2f]i)i/2(^g_^)i/2 ^^^Yi to = XQ/2ni in terms of the initial 
bubble size xq and thus exhibits a finite time singularity 
for X = 0, i.e., a zero bubble size or bubble closure at time 
to- In Fig. [2] we have depicted the finite-time-singularity 
solution for vanishing noise together with the noisy case. 



In the presence of thermal noise Eq. (|2.16p admits 
an exact solution, see e.g. Ref. [s^. The probabil- 
ity distribution, subject to the absorbing state condition 



III. WEAK NOISE ANALYSIS 

In the weak noise limit DkT — > we can apply a well- 
established canonical scheme to investigate the Fokker- 
Planck equation (|2.20p . see, for instance, Refs. p^. |20|. 
Introducing the WKB ansatz 



(3J) 



the weight (or action) S{x, t) satisfies the Hamilton- 
Jacobi equation 



dS_ 
'dt 



H ^0 



(3.2) 
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with Hamiltonian 



1 , / rii 



(3.3) 



From this scheme, the equations of motion yield in the 
form 



dx 
lit 
dp 



dt a;2 ^' 



(3.4) 
(3.5) 



They determine orbits in a canonical phase space 
spanned by the bubble size x and the momentum p. Com- 
paring the equation of motion (j3.4p with the Langcvin 
equation (|2.7p wc observe that the thermal noise ^ is re- 
placed by the momentum p = dS/dx. 

The action S associated with an orbit from xq to x 
during time t is given by 



S{x,t) 

or by insertion of Eq. (|3.4I 

S{x,t) = \ 



, dx 

dt p— - Ht, 

xo,0 



dt 



dt p^ 



(3.6) 



(3.7) 



xo,0 



A. Large bubbles 

For large bubbles, i.e., x ;» xi = ni/\Q2\, we can 
ignore the loop closure contribution characterized by fii, 
and we obtain the Hamiltonian 



H^^P^ + n^p, 

as well as the linear equations of motion 



dx 

'dt 
dp 
'di 



0. 



(3.8) 

(3.9) 
(3.10) 



The solution is given hy p — po^ x ~ xq + {po + )^ 
describing an orbit from (xo,Po) to {x,po) in time t. Iso- 
lating Po = {x — xo ~ ^2)/t and inserting in Eq. (|3.7p wc 
obtain the action 



1 {x - xq - n^tf 

^(^'*)=2 1 ' 



(3.11) 



and inserted in Eq. (j3.ip the biased random walk distri- 
bution ()2.1H) . In Fig.[3]wc have depicted the phase space 
for r^i — 0, i.e., in the large bubble-random walk case. 
The orbits are confined to the constant energy surfaces. 
We note in particular that the infinite time orbit lies on 
the p = —^2 manifold. We note, moreover, that in the 
large bubble case the weak noise case fortuitously yields 
the exact result for the distribution P. 
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Figure 3: We show the phase space structure in the case f2i = 
0, i.e., for random walk with constant drift. We show the 
zero energy manifolds for p = and p = —2Q2 and a negative 
energy orbit from xq to x in time t (arbitrary units). 



B. Small bubbles at and below Tm 

For small bubbles, i.e., x xi ~ ili/\^l2\, the loop 
closure contribution dominates and wc obtain the Hamil- 
tonian 



1 2 P^i 

2 X 



and the equations of motion 
dx 

'dt ~~ 
dp 

'dt~~ 



— +P, 

X 

—P, 



(3.12) 

(3.13) 
(3.14) 



determining orbits in {x,p) phase space. Eliminating p 
the bubble size is governed by the second order equation 



d X 
'd^ 

V = 



dV 

dx ' 

'2a;2' 



(3.15) 
(3.16) 



describing the 'fall to the center' {x = 0) of a bubble of 
size X, i.e., the absorbing state corresponding to bubble 
closure. 

The long time stochastic dynamics is here governed 
by the structure of the zero energy manifolds and fixed 
points. From Eq. (|3.12p it follows that the zero energy 
manifold has two branches: i) p = Q, corresponding to the 
noiseless transient behavior showing a finite time singu- 
larity as depicted in Fig. [2] and ii) p = 2ni/x associated 
with the noisy behavior. In Fig. [1] we have depicted the 
phase space structure. In the long time limit the orbit 
from xq to X passes close to the zero energy manifold 
p = 2f2i/x. Inserted in the equation of motion (|3.13p we 
have 



dx f2i 
dt x ' 
with long time solution 

x{t) - (217it)i/2. 



(3.17) 



(3.18) 
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p = 0, H=0 



the purpose of our discussion it is useful to introduce the 
parameters 



= c/2, 
7i / 1 



2k V T, 



(4.1a) 
(4.1b) 



Measuring time in units of /is the Fokker-Planck equation 
(|2.20p takes on the reduced form 

Note that ~ 1, and, close to the physiological temper- 
ature Tr, e « 2(T/T™ - 1). 



Figure 4: We show the phase space structure in the case ft2 = 
{T — Tm), i.e., for the small bubble dynamics governed by 
the entropic contribution. We show the zero energy manifolds 
p = and p = 2Qi/x and a negative energy orbit from xq to 
X in time t (arbitrary units). 



We notice that the motion on the noisy manifold p = 
2fli/x is time reversed of the motion on the noiseless 
manifold p = 0. Next inserting the zero energy manifold 
condition p = 2fli/x in Eq. (|3.7p we obtain 



S = 2ni / dt 



(3.19) 



and inserting the solution in Eq. p.lSp the action 

S{x,t) ^2Qilogx{t), (3.20) 
yielding according to Eq. (j3.1|) the long time distribution 
P{x,t) cx a;(17it)-^^i/2Z?fcT^ (3 2i) 

We have incorporated the absorbiiig state condition P = 
for a; = 0; as discussed in Ref. [3a| this condition fol- 
lows from carrying the WKB weak noise approximation 
to next asymptotic order. For the first-passage time den- 
sity of loop closure we obtain correspondingly 



W{t) (X t-''^/^DkT^ 



(3.22) 



We note that the power law dependence in Eqs. p.2ip 
and (|3?22)) is in accordance with Eqs. ((2T7)) and ((2?T8l) 
for DkT 0. 



IV. CASE OF ARBITRARY NOISE STRENGTH 

In the previous section we inferred weak noise-long 
time expressions for the distribution P on the basis of 
a canonical phase space approach. Here we address the 
Fokker-Planck equation (|2.20p in the general case. For 



A. Connection to the quantum Coulomb problem 

By means of the substitution P = e^-^x^''P, P satisfies 
the equation [2l[ 



dP _ 

'dt 2dx^ 



Id^P /^(Ai + l) 



2x2 



P, (4.3) 



which can be identified as an imaginary time Schrodinger 
equation for a particle with unit mass in the potential 



V{x)^ 



+ 1) _ , 
2x2 ^ 2 



(4.4) 



i.e., subject to the centrifugal barrier /i(/i + l)/x2 for an 
orbital state with angular momentum /i and a Coulomb 
potential — /ie/x. In Fig.[5]we have depicted the potential 
V — e2/2 in the two cases. In terms of the Hamiltonian 



H 



1 <P /i(^ +1) lie £2 

2rfx2 2.x2 X 2 



(4.5) 



the eigenvalue associated with Eq. (|4.3p problem has the 
form 



(4.6) 



Expressed in terms of the eigenfunctions the transition 
probability P(x,xo,t) then becomes 

P(x,xo,t) = e^"-^") (^)''^e-^"**„(x)vI/„(xo). 

n 

(4.7) 

Here, the completeness of ensures the initial condition 
P(x, xo, 0) = 5{x — xq). Moreover, in order to account for 
the absorbing boundary condition for vanishing bubble 
size we choose \l/n(0) = 0. We also note that for a finite 
strand of length L, i.e., a maximum bubble size of L, we 
have in addition the absorbing condition 4'„(L) = for 
complete denaturation. Expression (|4.7p is the basis for 
our discussion of DNA-breathing, relating the dynamics 
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the bubble dynamics, corresponding to complete denat- 
uration of the DNA chain. 

Mathematically, we model the bubble dynamics with 
absorbing boundary conditions at zero bubble size x = Q, 
and, for a finite chain of length L, a± x = L. When 
the bubble vanishes or complete denaturation is reached, 
that is, the dynamics stops. Physically, this stems from 
the observation that on complete annihilation (closure) 
of the bubble, the large bubble initiation barrier pre- 
vents immediate reopening of the bubble. Similarly, a 
completely denatured DNA needs to re-establish bonds 
between bases, a comparatively slow diffusion-reaction 
process. 



1. Long times for T < Tm 

At long times and fixed x and xq, it follows from Eq. 
(|4.7p that the transition probability is controlled by the 
bottom of the energy spectrum. Below and at Tm the 
spectrum is continuous with lower bound Setting 
Ek = -I- in terms of the wavenumber k and 
noting from the eigenvalue problem in Eqs. (|4.5p and 
(1461) that -^k^x) ^ for small kx we find 

P{x,XQ,t) oc cxp {~\e\{x - xo)) (^^^ (^~2") 

(4.8) 



Figure 5: Schematic of the potential V{x) — /2. a) T < Tm'. 
The potential is repulsive, yielding a continuous spectrum. 
The bubble fluctuations correspond to a biased Brownian 
walk process in bubble size x before collapse at a; = 0. b) 
T < Tm- The potential is attractive and can trap a series of 
bound states. At long times the lowest bound state indicated 
in the figure controls the behavior. The bubbles increase in 
size eventually leading to complete denaturation. 



to the spectrum of eigenstates, i.e., the bound and scat- 
tering states of the corresponding Coulomb problem [s^ . 

The transition probability P{x,XQ,t) for the occur- 
rence of a DNA bubble of size x at time t is controlled by 
the Coulomb spectrum. Below the melting temperature 
Tm for e cx {T/Tm — 1) < 0, the Coulomb problem is re- 
pulsive and the states form a continuum, corresponding 
to a random walk in bubble size terminating in bubble 
closure {x = 0). At the melting temperature Tm for 
e = 0, the Coulomb potential is absent and the con- 
tinuum of states is governed by the centrifugal barrier 
alone, including the limiting case of a regular random 
walk. Above the melting temperature for e > 0, the 
Coulomb potential is attractive and can trap an infinity 
of bound states; at long times it follows from Eq. (|4.7p 
that the lowest bound state in the spectrum dominates 



By a simple scaling argument we then obtain the long 
time expression for the probability distribution 

P{x, xo,t) (X x4+''^e-l^l(^-"°)e-^'*/2f-3/2-M. (4.9) 

The lifetime of a bubble of initial size xq created at time 
t = follows from Eq. (|4.9p by calculating the first pass- 
ing time density W{t) in Eq. (|2.13p . Using the Fokker- 
Planck equation (|4.2p we also have more conveniently 



|^+(?^-2.1P 

ox \ X 



and we obtain at long times 



W{t) cx (1 + 2/^)a;^+^''el^l^«e 



\e\x(.-e^t/2,-3/2-ti 



(4.10) 



(4.11) 



In Fig. [HI we have depicted the bubble lifetime distri- 
bution W{t) below Tm for e = -1/2. 



2. At the transition T ^ Tm (e ^ 0) 

At the transition temperature T ~ r,„ for e = the 
Coulomb term is absent and we have a free particle sub- 
ject to the centrifugal barrier //(/i -|- l)/2x^. In this case 
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Figure 6: Bubble lifetime distribution W{t) from Eq. (|4.11|l . 
with e = -1/2, xo = 5, and c = 1.76 (full line) and 2.12 
(dashed). The initial power-law behavior with slopes -2.38 
and -2.56 is indicated by the straight lines. Inset: log versus 
linear scale, emphasizing the exponential decay for long times. 



t igure 7: Bubble lifetime distribution W{t) from Eq. (|4?T6)) 
for T = T,n, Xo ~ 5, as well as c = 1.76 (full line) and c = 2.12 
(dashed line). Inset: log-log plot of the power-law behavior 
at long t, with slopes —2.38 and —2.56, as indicated by the 
straight lines. 



the eigenfunctions are given by the Bessel function [36[ 
^kix) = {kxf/^Ji,2+^.{kx), (4.12a) 



Ek = — , 



(4.12b) 



where orthogonahty and completeness follow from the 
Fourier-Bessel integral [s^ 



kJ,{kx)dk / yj,{ky)f{y)dv (4.13) 



By insertion into Eq. (j4.7p we obtain the distribution 
P{x,xo,t) 



j'OO 

(fcx)Ji/2+^(fcxo), 

"'0 



or, by means of the identity [3f 



Uax)Uhx)xdx = (^g 



the explicit expression 



(4.14) 



P{x,x,,t) = (^)"(x:ro)i/2f~ie-(-^+-o)/2* 



'xh/2+tL{xxo/t). 



(4.15) 



Here, Iv{z) is the Bessel function of imaginary argument 
[set . From Eq. (|4.15p we also infer, using Eq. (|4.10p the 
first passage time distribution 



W{t) 



2x. 



r(i/2 + M) 



^o/^*(2t)"^/^"^. 



(4.16) 



in accordance with Eq. (|4.1ip for e = 0. In Fig.[7]we show 
the first passage time distribution (|4.16p for two different 
critical exponents c. Note that the power-law exponent 
—3/2 — u = — 3/2 — c/2 is identical to the result reported 
in Ref. 



3. Long times for T > Tm 

Above the transition temperature for e > the 
Coulomb potential — /le/x is attractive and can trap a 
series of bound states. In the long time limit the low- 
est bound state controls the behavior of P. According 
to Eqs. (|4.5p and (|4.6p the lowest bound state with 
eigenvalue Ei < e^/2 must satisfy the eigenvalue equa- 
tion 



Id^^i /i(/i + l) 

2 d.T2 ^ 2x2 ~ 



*i = Si*!. (4.17) 



For a; -> cx) we have -(1/2)*'/ = {Ei - eV2)^i 
and 'I'l must fall off exponentially, 'I'l ~ cxp(— Ax), 
A = (2£;i - e2)i/2. For X ^ we have -(1/2)*'/ + 
(/i(ju + l)/2x2)vl/i ~ and we infer *i ^ x^+^. Conse- 
quently, searching for a nodeless bound state of the form 
*i ^ x^+'^ exp(— Ax) we readily obtain the normalized 
lowest level 

*i(x) = ylxi+'^e-'^'^/(i+'^\ (4.18a) 
(2/ie/(/i+l))2A'+3 

^ - r(2, + 3) ' (^-^^^^ 

with corresponding eigenvalue 

i?i = ^(l-(M/(A* + l))'). (4.19) 
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The maximum of the bound state is located at (/i + 
l)^/^e ~ l/(r — Tjn) and thus recedes to infinity as we 
approach the melting temperature. From Eq. (|4.7p we 
thus obtain after some reduction 

^g-e^(l+2M)t/2(l+M)2^ (4.20) 

Above T„i the bubble size, on average, increases in time 
until full dcnaturation is reached. In terms of the free 
energy plot in Fig. [T|d this corresponds to a Kramers es- 
cape across the (soft) potential barrier (corresponding to 
a nucleation process). This implies that the transition 
probability P{x,xo,t) from an initial bubble size xq to 
a final bubble size x must vanish in the limit of large t. 
According to Eq. (|4.20p P{x,xo,t) decays exponentially, 



P{x, xo,t) (X e 



with a time constant given by 
2(1 + m)^ 



cx IT - T„ 



(4.21) 



(4.22) 



Exact results 



The eigenvalue problem given by Eqs. (|4.5p and (|4.6 



1 (P-^ 
2 



1) _ M£ 
2a;2 X 



* = (4.23) 



has the same form as the differential equation satisfied 
by the Whittaker function w [3^, 



dx^ 



A 1/4 -m^ 



w = 0, 



(4.24) 



with the identifications z = 2kx, A — fJ^e/n, m = 1/2 + fi^ 
and E = e'^/2 — /2. Incorporating the absorbing state 
condition \E'(0) = and using an integral representation 
for the Whittaker function w (Sal we obtain the solution 



'^{x) cx (2ra)^+^e-''^ 

(4.25) 



X / e 
/o 



In the bound state case for e > the parameter k > 
and the bound state spectrum is obtained by terminating 
the power series expansion of Eq. (|4.25p 38| , 



^{x) cx (2Ata;)^+^e-''^ 

x$(l + /i(l - e/k), 2(1 + ^); 2kx), (4.26) 

with the polynomial 

a z a{a. + 1) z^ 



$(a,7;z) = 1 



7I! 7(7 + 1)2! 
a(a + l)(a + 2) z^ 
'7(7 + l)(7 + 2) 3!- 



(4.27) 



Simple algebra then yields the spectrum 
K = e— y— , n = l,2,--- 

and associated eigenfunctions 

^ cx x^^^e^'^^ X polynomial, 



(4.28) 



(4.29) 



the lowest state and eigenfunctions given by Eqs. (|4.18ap 
and (jil^ . 



V. DISCUSSION 



In typical experiments measuring fluorescence corre- 
lations of a tagged base pair bubble breathing can be 
measured on the level of a single DNA molecule . 
The correlation function C (<) is proportional to the inte- 
grated survival probability, i.e., 



C(t) cx / P{x,xo,t)dx, 
Jo 



(5.1) 



where L is the chain length. From the definition of the 
first passage time distribution in Eq. (|2.13|) we also have 



C{t) = 1 - / Wit')dt'. 
Jo 



A. Below T„, for e < 



(5.2) 



Below the melting temperature T!,„ < wc obtain from 
Eq. gHI) 

C(t) = 1 - 4+'^el^l"° f e-'"''/\t')-^/^~''dt\ (5.3) 



or in terms of the incomplete Gamma function 7(0;, x) = 
e-~H°"^dt [Hi 

C{t) = l-4^^''el'l^"(eV2)i/2+M 

X7(-l/2-Ai,eV2). (5.4) 

Using 7(a, x) — T{a) — x"~^e~^ for x 00 we have for 
large t 

Cit) = const. + 4+2'^e-2el^l-°f-3/2-Me--^t/2. (5.5) 

We note that the basic time scale of the correlations is 
set by cx (r,„ — T)^^. As we approach the time 
scale diverges like (T^ — T)~^. 

For t ^ the correlations show a power law behavior 

Cit) = const. + C(t) cx i"^/2~''(mod a const.), (5.6) 

with scaling exponent —3/2 — /i = —3/2 — c/2. Here 
3/2 originates from unbiased bubble size random walk 
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Figure 8: The structure function C{ijj)- For |ci;j <c; ^ the 
structure function has a Lorentzian Une shape; for ^ 
it exhibits power law tails. 



whereas the contribution /.j = c/2 is associated with the 
entropy loss of a closed polymer loop. 

At long times t ^ the correlations fall off expo- 
nentially 

C{t) = const. + C{t) oc e"''*/^(mod a const.). (5.7) 

The size of the time window showing power law behav- 
ior increases as Tm is approached. This corresponds to 
the critical slowing down on denaturation, as already ob- 
served in Ref. [lj| numerically, and in Rcf. [l7[ in absence 
of the critical exponent c due to polymeric interactions. 
In frequency space the structure function is given by 



^C{t)dt. 



(5.8) 



By means of a simple scaling argument we infer that C{uj) 
has a Lorentzian line shape for |a;| <C crossing over to 
power law tails for 3> e^. 



C{u) - Xq 



l+2M„|e|a;o 



tj2 + (e2/2)2 



1 



1/2+M 



for |cj| < e2(5.9a) 
for > e=^(5.9b) 



In Fig. [5] we have depicted the structure function C{uj). 



B. At for e = 

At the transition temperature Tm the exact expression 
for the first passage time distribution is given by Eq. 
dHH). Using Eq. for C{t) we then obtain 



C[t) 



1 



T{l/2 + ii,xl/2t) 
r(l/2 + M) ' 



(5.10) 



where T{a^) = e *t" ^dt is the incomplete Gamma 
function 13611. 



At short times we have 
C{t) = 1 
whereas for t 

C{t)^ 



r(i/2 + M) 

2(xg)^/2+^' 



2t 



^1/2-/^5-2:^/ 



^-l/2-/._ 



(5.11) 



(5.12) 



The correlation function thus exhibits a power law be- 
havior with scaling exponent —1/2 — /i = — 1/2 — c/2, as 
obtained from a different argument in Ref. [29|. Corre- 
spondingly, the structure function C{lj) has the form 



(5.13) 



(7(u;)cx4+'^|a;P'-i/2. 
C. Above T,„ for e > 



Above Tm (e > 0) the DNA chain eventually fully de- 
natures and the correlations diverge in the thermody- 
namic limit. We can, however, at long times estimate 
the size dependence for a chain of length L. From the 
general expression (|4.7p we find 



C{t) 







(5.14) 



At long times the lowest bound state dominates the ex- 
pression. Inserting ^'i and E'l fromEqs. (j4.18ap . (|4.18bp . 
and (j4.19p and performing the integration over x we ob- 
tain 



C{t) oc A^e— 

X (1 + ii)€-^ [1 + {Le/{\ + l)e^^/(i+^)]] . (5.15) 

The correlations decay exponentially with time constant 
~ e~'^{lj,+ l)2/(2/i-|- 1). In frequency space the structure 
function has a Lorentzian lineshape of width e^(2/x -|- 
l)/(/i + 1)2, and for the size dependence one obtains 



C{t) 



^geL/(i+,0^ for eL/(l + ^) > 1, 
Le/{1 + fi), for eL/{l -|- ^) < 1 



(5.16) 



Note that close to Tm the correlation function C{t) cx 
L. In Fig. [H] we depict in a plot of C/L vs. L the size 
dependence of the correlation function. 



D. Comparison to experimental data 

Below the melting temperature Tm, DNA breathing 
can be monitored on the single DNA level by fluores- 
cence correlation spectroscopy 0, [13, [l^. In the FCS 
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C/L 



exponential 



constant 



{^+^J)/E 



Figure 9; We depict C/L as a function of L. For L ^ {l+fi)/e 
the correlations depends linearly on L; for L 3> (1 + /i)/e the 
correlations increase exponentially as a function of L. 



FCS data (AT)„, T=55°C 
FCS data A, o, T=20°C 
FCS dataM,,, T=30°C 



r model, ^=0 
r model, ^=0, 
r model, ^=1 , 



=2""\x„=0.692 



0.24 

3.6""^, Xo=1.08 
Diffusion model 
Diffusion model, long-time limit 



"■S 0.2 



o 0.1 
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Figure 10: (Color online). Drift-diffusion model and experi- 
mental data from Ref.[^ compared to the F model, for vari- 
ous parameters. The curve for jj, = and e = l/\/2 exactly 
matches the long-time behavior from Ref. 



experiment from Ref. a DNA construct of the form 



5' I T 

GGCGCCCATATATATATATATATATATGCGC T 
CCGCGGGTATATATATATATATATATACGCG T ' ' 

3' I T 

Q 

was employed. Here, a bubble domain consisting of 
weaker AT base pairs are clamped by stronger GC base 
pairs. On the right, a short loop consisting of four 
T nucleotides is introduced. The fluorophore (F) and 
quencher (Q) are attached to T nucleotides as shown. 
With the highest probability, a bubble will form in the 
AT-bubble domain. As the bubbles consist of flexible 
single-strand, in an open bubble the fluorophore and 
quencher move away from each other, and fluorescence 



occurs. Once in the focal volume of the FCS setup, bub- 
ble opening and closing corresponds to blinking events 
in the signal, whose correlation function (corrected for 
the diffusion in and out of the focal volume) are shown 
in Fig. [TOl Three different bubble domains with chang- 
ing sequence were used to check that potential secondary 
structure formation does not influence the breathing dy- 
namics, confirming the picture of base pair-after-base 
pair zipping and unzipping. The figure shows examples 
from all three constructs, underlining the data collapse 
already observed in Ref. Q. 

The theoretical lines shown in Fig. [TO] correspond to 
the biased diffusion model introduced in the original ar- 
ticle [1]. While the full solution of this diffusion model 
fits the data well over the entire window, the long time 
expansion demonstrates the rather weak convergence of 
the expansion. In Fig. [TO] we also included our asymp- 
totic solution (|5.4p for the autocorrelation function, for 
various parameters. Good agreement with the data is 
observed. 



VI. SUMMARY AND CONCLUSION 

In this paper we have analyzed the breathing dynam- 
ics of thermally induced denaturation bubbles forming 
spontaneously in double-stranded DNA. We have shown 
that the Fokker-Planck equation can be analyzed from 
two points of view: i) In the weak noise or low temper- 
ature limit a canonical phase space approach interprets 
the stochastic dynamics in terms of a deterministic 'clas- 
sical' picture and gives by simple estimates access to the 
long time dynamics. In particular, we deduce that the 
dynamics at the transition temperature is characterized 
by power law behavior with scaling exponent depend- 
ing on the entropic term, ii) In the general case we show 
that the Fokker-Planck equation can be mapped onto the 
imaginary time Schrodinger equation for a particle in a 
Coulomb potential. The low temperature region below 
the transition temperature then corresponds to the con- 
tinuum states of a repulsive Coulomb potential, whereas 
the region above Tm is controlled by the lowest bound 
state in an attractive Coulomb potential. The mapping, 
moreover, allows us to calculate the distribution of bub- 
ble lifetimes and the associated correlation functions, be- 
low, at, and above the melting temperature of the DNA 
helix-coil transition. Finally, at the melting transition, 
the DNA bubble-breathing was revealed to correspond 
to a one-dimensional finite time singularity. 

The analysis reveals non-trivial scaling of the first pas- 
sage time density quantifying the survival of a bubble 
after its original nucleation. The associated critical expo- 
nent depends on the parameter ^ = c/2 stemming from 
the entropy loss factor of the flexible bubble. The first 
passage time distribution and correlations depend on the 
difference T/Tm — 1, and therefore explicitly on the melt- 
ing temperature Tm (and thus the relative content of AT 
or GC base pairs). We also obtained the critical depen- 
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dence of the characteristic time scales of bubble survival 
and correlations on the difference T—T^. The finite size- 
dependence of the correlation function was recovered, as 
well. 

The mapping of the of DNA-breathing onto the quan- 
tum Coulomb problem provides a new way to investigate 
its physical properties, in particular, in the range above 
the melting transition, T > T^. The detailed study of 
the DNA bubble breathing problem is of particular inter- 
est as the bubble dynamics provides a test case for new 
approaches in small scale statistical mechanical systems 
where the fluctuations of DNA bubbles are accessible on 
the single molecule level in real time. 
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